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Abstract 

Ab-initio full-potential linearized augmented plane wave (FLAPW) calcula- 
tions have been used to study the influence of the interface morphology and, 
notably, of the exchange reaction on the electronic properties of Al/GaN (100) 
interfaces. Although the detailed mechanism is not understood, the exchange 
reaction has been purported to influence the Schottky barrier height as a result 
of the formation of an interfacial Ga^^Ali-^^N layer. In particular, the effects 
of interface structure (i.e. interfacial bond lengths, semiconductor surface po- 



larity, and reacted intralayers) on the SBH at the Al/GaN (001) junction are 
specificahy addressed. Thus, the electronic structure of the following atomic 
configurations have been investigated theoretically: (i) an abrupt, relaxed 
GaN/Al interface; (ii) an interface which has undergone one monolayer of ex- 
change reaction; and interfaces with a monolayer-thick interlayer of (in) AIN 
and (iv) Gao.5Alo.5N. The exchange reaction is found to be exothermic with 
an enthalpy of 0.1 eV/atom. We find that the first few layers of semicon- 
ductor are metallic due to the tailing of metal-induced- gap-states (MIGS); 
therefore, the presence of a monolayer-thick interfacial alloy layer does not 
result in an enhanced bandgap near the interface. Intermixed interfaces are 
found to pin the interface Fermi level at a position not significantly different 
from that of an abrupt interface. Our calculations also show that the interface 
band line-up is not strongly dependent on the interface morphology changes 
studied. The p-type SBH is reduced by less than 0.1 eV if the GaN surface is 
Ga-terminated compared to the N-terminated one. Moreover, we show that 
both an ultrathin Gax.Ali_i.N (x = 0, 0.5) intralayer and a Ga<-^A1 atomic 
swap at the interface does not significantly affect the Schottky barrier height. 
PACS 71, 73, 73.40.Sx, 73.61.-r 



I. INTRODUCTION 

In the last few decades, a significant number of studies have investigated 
metal/semiconductor interfaces, due to the key role played by ohmic and Schottky con- 
tacts in technological III-V semiconductor devices. However, from both the experimental 
and the theoretical points of view, the fundamental mechanism involved in the Schottky 
barrier formation has not been fully elucidated^Fi. 

Al / III-V semiconductor junctions have been the most extensively studied interfaces to 
date. This has occurred due to both practical and fundamental reasons. Al contacts are 
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inexpensive to manufacture and have stable electrical, chemical and structural characteristics 
at moderate temperaturesl. Al contacts also are of practical interest since they are commonly 
used in IIl-V commercial device technology, particularly as the gate of field effect devices. 

In the field of ab-initio computational physics, many works have addressed band align- 
ment at semiconductor heterojunctionsi, whereas few calculations have studied the barrier 
height of metal/semiconductor interfaces. Al has been used almost exclusively as the metal 
layer in theoretical investigations of metal/semiconductor contacts0^0. This arose because 
Al is found to be advantageous since it lacks d-electrons, forms a free-electron metal, and 
lattice matches a number of Ill-V semiconductors. Recently, the GW method has been used 
to study the Al/ GaAs interfaceS. 

More recently, theoretical work on III-V interfaces demonstrated that the pressure- 
induced barrier height changes can be used to critically determine the nature of the 
states which pin metal/III-V interface^llli. Experimental measurements of the pressure- 
induced change in Al/GaAs and Al/AlGaAs barrier heights are consistent with defect- free 
interfaces0'0, validating Al/III-V interfaces as a model system. In contrast, Au, a metal 
which reacts with GaAs to release near-interfacial As at GaAs interfaces, exhibits pressure- 
induced barrier height changes consistent with an interface decorated with a deep-level 
point-defect, such as the As antisite (Asca)^- Because GaN solid is significantly more sta- 
ble than GaAs, the Fermi- level pinning position of Al/GaN interfaces would also be expected 
to be determined by the properties of interfaces without deep-level defects. 

In the study reported here, we start an investigation of an important aspect of Schottky 
barriers that only recently has begun to be explored using modern electronic structure 
calculations: the role that interface chemical exchange plays in the electrical and electronic 
properties of interfaces. Al/GaN interfaces were used in this study due to the well known 
chemistry and atomic structure of the interface. GaN0, GaAs0, and InPta, among others, 
are found to exhibit an exchange reaction in which Al replaces the Group III element at 
the surface. The extent of the reaction is found to be of the order of a monolayer for 
unannealed contacts and increases at elevated annealing temperatureJii. Changes in the 



electrical properties, as reported in a large number of III-V studies, have attributed the 
increase in the barrier for n-type semiconductors to the increased bandga p@. However, a 
more recent study on both n-type and p-type contacts indicates that the observed behavior 
can be attributed to a shift in interfacial Fermi level pinning position, rather than due to 
the formation of an increase in interfacial bandg apil. A strong fundamental understanding 
of the influence of the chemical reaction on the interface electronic structure has not yet 
been established. 

Density functional calculations are able to quantitatively address many of the ques- 
tions including the energetics of the exchange reaction, the nature and extent of the MIGS 
(metal-induced gap states) and the role of the larger-bandgap AlGaN interface layer on 
the interface electronic structure. In a previous workil, we performed ab-initio calculations 
to determine the electronic properties of the [001] ordered atomically abrupt N-terminated 
XN/Al (X = Ga, Al) interface, focusing mainly on the Schottky Barrier Height and the 
resulting interface states. In the work reported here, the structural and electronic properties 
of several Al/GaN interface configurations are investigated. In particular, the effects of in- 
terface structure {i.e. interfacial bond lengths, semiconductor surface polarity, and reacted 
intralayers) on the SBH at the Al/GaN (001) junction are specifically addressed. It should 
be noted that some experimental results^ are available for the (0001) wurtzite GaN/Al in- 
terface and that some caution has to be taken in comparing our results with the experiments, 
since the polarization effects and the different coordinations of the surface bonds may both 
play a role in determining the structural and electronic properties. 

Starting from the atomically abrupt N-terminated GaN/Al junction (in the following 
denoted as the A structure), we consider four other different interface morphologies: (i) the 
ideal, abrupt, relaxed Ga- terminated GaN/Al interface (B); (ii) the configuration with an 
AIN intralayer (C); (Hi) the configuration in which an Al atom has undergone an exchange 
reaction with a surface Ga atom (D) and finally (iv) the configuration with a Gao.5Alo.5N 
alloy layer at the interface (E). Our calculations show that the interface band line-up is not 
strongly dependent on the interface morphology. In particular, we find that the exchange 
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reaction is found to be exothermic with an enthalpy of 0.1 eV/atom. We show that both 
an ultrathin Gaa.Ali„j.N (x= 0, 0.5) intralayer and a Ga^Al atomic swap at the interface 
does not significantly affect the p-type SBH. The p-type SBH is smaller by 0.1 eV for a 
Ga-terminated surface than a N-terminated one. 

In Sec. |T|, we describe the theoretical method and the atomic configurations used in 



this study. In Sec. |ITl[ we discuss our results, focusing in Subsec. [Ill A| on the structural 
properties, {i.e. equilibrium atomic distances, nature of the interface bonds - metallic 
versus covalent) whereas in Subsec. [1IIB| we discuss effects of the interface morphology on 



the interface states and the resulting SBH. Finally, our results are summarized in Sec. 1^. 

II. COMPUTATIONAL METHOD AND ATOMIC CONFIGURATIONS USED IN 

CALCULATIONS 

First principles full-potential linearized augmented plane wave (FLAPwi) calculations 
within the local density approximation to density functional theory were used in this work@. 
Plane waves with wave vector up to Kmax = 4.6 a.u., leading to about 5000 basis functions, 
were used. Angular momenta up to Imax = 8 in the muffin tin spheres (-Rca =Rai = 2.1 a.u. 
i?N = 1-4 a.u.) for both wavefunctions and charge density were used in the self-consistency 
cycles. Three special k-points, obtained following the Monkhorst-Pack schemeil, were used 
to sample the tetragonal Brillouin zone. 

The widely adopted supercell approach was used to simulate the different types of inter- 
faces. In particular, we performed ab-initio calculations for supercells containing 15 GaN + 
7 Al (13 GaN + 7 Al) layers in the N-terminated (Ga-terminated) case. The thickness of 
these layers is sufficiently large that the metal and semiconductor layers farthest from the 
interfaces have properties similar to the bulk (see below). 

The method employed to evaluate the SBH uses the metal and semiconductor atomic 
core levels as reference energies located sufficiently far from the interface to have bulk- 
like bandstructure. In particular, the SBH can be expressed as $n = Ab + AEh, where 
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Ab and AEi, indicate an "interface" and "bulk" contribution, respectively (see Ref. ^ for 
details in the GaN/Al case). This simple procedure is commonly used in ab-initio all- 
electron calculations to evaluate the interface band line-up at semiconductor heterojunc- 
tions or metal/semiconductor junctions. We considered all the structures for a cubic {i.e. 
zincblende) [001] ordered GaN substrate, with a calculated lattice constant asubs = ciCaN = 



8.47 a.u.E5l. Further structural details can be found in Ref. ^ All the structures considered 
are shown in Fig. |l] and summarized in Table |l[ 

The choice of these structures (in particular the C, D and E systems) has been sug- 
gested by a recent surface investigation of Al films deposited on wurtzite GaN(0001)-(lxl) 
surfaces,@ reporting that the interface is very reactive. It was experimentally observed, in 
fact, that metallic-Ga was released from the interface after monolayer depositions of Al. 
This suggests that an exchange reaction occurs in which the Al atom replaces the group III 
element within the semiconductor, in analogy to that found for GaAs0 and InPEl. After 
repeated cycles of Al deposition and annealing, evidence was found0 for the presence of a 
Ga^jAli.^N alloy interfacial layer. 

indicated by 

In Table | (structure E) the interface plane X, indicates the fee cationic sites are occupied 
either by Ga or Al (50 % of Ga and 50 % of Al), with an ordered superstructure which aims to 
simulate the Gao.5Alo.5N alloy. In this case, we consider a supercell with doubled dimensions 
in the x — y plane. The resulting increase in the computational cost was compensated for 
by reducing the thickness of the supercell {i.e. 11 GaN + 5 Al layers). This is expected to 
affect the final value of the SBH, since bulk conditions in the thinner GaN and Al regions 
may not be completely recovered (see below). However, the effects of the alloy Gao.5Alo.5N 
intralayer on the structural and electronic properties can be determined by comparing the 
results to the ideal atomically abrupt N-terminated GaN/Al interface having the same size 
and in-plane dimensions; the only difference between them being the last cationic plane of 
the semiconductor side which is made of 50% Ga and 50% Al atoms and all Ga atoms, in 
the first and second case, respectively. 



For configurations A-D, structural relaxations along the [001] direction were performed 
according to ab-initio atomic forcesH. Other degrees of freedom including in-plane atomic 
relaxations and formation of extended defects, such as dislocations, were neglected. In 
the case of configuration E, due to onerous computational costs for the in-plane doubled 
supercells, structural minimization was not performed. Instead, the calculated interplanar 
distance from structure A was used. Despite this limitation, the comparison of the results 
for configurations A-E can be used to determine the dependence of the interfacial electronic 
properties on the chemistry and structure of the interface. 



III. RESULTS AND DISCUSSION 



A. Stability and structural properties 



with the supercells 

Since the supercells do not have the same number of each atomic species, a fully consistent 
determination of the relative stability of the configurations is in general not possible from 
the superlattice calculations alone@. The problem can be avoided by obtaining the chemical 
potentials /i of the atomic species of pure solids from independent total energy calculations. 
However, the problem of biggest interest in the present contest is the one of the exchange 
reaction, and we can get the relevant figure by comparing the ideal N-terminated interface 
(A) with the structure showing the Ga<-^A1 swap (D), since the two supercells contain the 
same number of atoms of the same species. As for the Ga vs N termination, this problem 
is of more concern in the study of the GaN surface, and has already been addressed in the 
literatureii. As shown in Ref. and mentioned earlier in this work, the atomic forces acting 
on the interface Alj in configuration A do not decrease the zai,-n^ pushing the Al towards 
the GaN region. Instead, we find an equilibrium structure having a bond length dA^-Ni 
larger than both d^Xi-N do'l^'tj^ (where d'^x-N 0^ = Ga, Al) denotes the bulk bond length 
of XN). However, our calculations show that the total energy of configuration D is lower 
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(by about 0.1 eV/atom) than the total energy of structure A, confirming that the Ga<-s>Al 
exchange reaction observed experimentallyEl is energetically favorable. Intuitively, this can 
be understood by noting that the AIN bond has a significantly higher bond energy than 
GaN. We note that a further displacement of Ga inside the Al overlayer will imply a more 
ionic N-Al interface bond (greater cohesive energy), and one more partially covalent Ga-Al 
bond on the metal side. We can therefore reasonably expect that such an exchange will 
lead to a further total energy reduction. The heat of formation for the Al/GaN exchange 
reaction is significantly less than the values predicted for the Al/GaAs system by Ihm (0.48 
eV)§ and Swarts (0.62 eV)0. 

The theoretically-determined interplanar distances along the [001] direction are summur- 
ized in Table |I[ Since the electronic properties of configuration A have been discussed in 
a previous workEl, we describe these results only in comparison with the other configura- 
tions. In the atomically abrupt interface B, the interplanar distance in the last GaN plane 
{/\zGai-Ni = 2.23 a.u.) is larger than in bulk GaN (2.12 a.u.), which indicates that the 
interface GajNj covalent bond is weakened with respect to the GaN bond in the bulk. This 
can be explained in terms of the half-metal/half-semiconductor character of the interface 
Gaj atom. In fact, as already pointed out in Ref. ^ for the A configuration, the interface 
cation {i.e Alj and Ga^ in the A and B case, respectively) forms with the other Al atoms 
a metallic bond, which reduces its possibility of full sp^ hybridization, thus weakening the 
covalent character of the N-cation bond. A similar mechanism is also responsible for making 
the interplanar distance Az;v2^;2 larger for the C configuration than Azj<^iij^2. In order to 

i i i i 

substantiate this, we report in Fig. |^ (a) and (b) the charge density contour plots relative 
to the B configuration for the bulk Al atoms and the interface Gaj-Alj atoms, respectively. 
The interface Gaj-Alj bond is predominantly metallic, although it is slightly different from 
the pure Al-Al bond, due to different Ga and Al electronegativities (compare Fig. (a) and 
(b)). 

In Fig.^ (c) we report the charge density contour plots for the interface Gaf and Al^ atoms 
in configuration D; the similarity of panels (b) and (c) in Fig. 1 shows that the Ga atom 



released in the swap mechanism tends to forms metalhc bond with the deposited Al atoms. 
This is consistent with the resuhs from ab-initio atomic forces, which give an interplanar 
distance for interface D nearly equal (within 1%) to that in the B configuration (see Table 
|). We note that the results obtained in the present work for configuration B are similar to 
those obtained from ab-initio pseudopotential calculations performed for the Ga-terminated 
Al/GaAs(001) interfaced the top Ga atomic layer was found to relax outwards, thus 
elongating the interface Ga-As bonds, as a result of metallization of the Ga-Al layer. 

Table I shows that in the fully relaxed C configuration, the Nj - Alj and Alj - Nf inter- 
planar distances (A2;^;i_^. ~ 2.12 a.u.) are equal to the interplanar distance {Azcai-Ni ~ 
2.12 a.u.) in configuration A. On the other hand, the LDA-calculated equilibrium interpla- 
nar distance for tetragonal AIN with the in-plane lattice constant fixed to that of the GaN 
substrate is Az^^i^"'^ = 2.00 a.u., leading to a bond length d^Ai^N = 3.60 a.u..0 

Therefore, our structural results regarding configuration C suggest that the AIN in- 
tralayer, constrained on a GaN substrate, can not attain the equilibrium tetragonal bond 
length within one layer. Instead, when Al replaces the surface Ga atom, the Al is found to 
occupy the same position without significant structural relaxation. Based on these results, 
we can reasonably expect a similar situation in the case of configuration E, with the Al 
atoms of the Gao.5Alo.5N layer substituting 50 % of the Ga atomic sites without appreciably 
modifying the bond lengths, thus excluding the possibility of buckling effects. Hence, the 
unrelaxed E configuration, whose interplanar distances are set to those of configuration A, 
should not be radically different from the fully relaxed structure. 

B. Schottky Barrier Heights 

This subsection focuses on the influence of the interface morphology on the SBH. We 
recall that for the GaN/Al interface we found that the metal induced gap states (MIGS) tail 
into the semiconductor side with a decay length of about A ~ 3.5 a.u.ii; in fact, in the present 
cases we flnd that the flrst few layers of the semiconductor have metallic character. For this 
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reason, the presence of a monolayer-thick interfacial alloy layer does not result in an enhanced 
bandgap near the interface. As a consequence, in this case changes in the electrical properties 
arise from shifts in the interface Fermi-level pinning position within the GaN bandgap. 
Table |I| lists our theoretical results for the p-type SBH = Ep — Evbm calculated as 
the difference between the Fermi level (Ep) and the GaN valence band maximum (Evbm) 
sufficiently far from the interface. Note that the n-type SBH, can be obtained from 
the p-type SBH through the following relation 

Egap being the experimental semiconductor band gap energy (E^g^piGaN) = 3.39 evi'0). 

We point out that all the results reported in the first row of Table |I| refer to supercells 
with 15 GaN + 7 Al layers. In order to check the convergence of our SBH values as a function 
of the cell dimensions, we also performed calculations for a larger supercell, containing 19 
GaN + 9 Al layers. /^From this caseS, we obtained variations of the SBH of less than 0.02 
eV, showing that bulk conditions are well recovered also in the smaller 15 GaN + 7 Al 
supercell. Note also that for configuration A we used an even smaller cell size (11 GaN + 5 
Al layers - see the second row in Table ^ in order to compare with configuration E, which, 
due to the in-plane doubling of the cell, was calculated with a shorter periodicity along the 
growth direction. The different value (by about 0.2 eV) obtained for the 15 GaN + 7 Al 
(configuration A, first row) with respect to the 11 GaN + 5 Al (configuration A, second row) 
is related to the insufficient thicknesses considered. However, this is not going to affect our 
conclusions, since in this case we are interested in SBH trends rather than absolute values. 
The results shown in Table |I| do not include quasi-particle corrections (Aqp ^0.1 eVi) 
related to the difference in screening on the two sides of the interface. This correction should 
be added to the p-type SBH value before a direct comparison to experiment is made. 

Our calculations show that the interface band line-up is not strongly dependent on 
the interface morphology and that all the SBH values found are in qualitative agreement 
with the experimental value reported for the atomically clean wurtzite-GaN/Al interface^. 
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We first discuss the results obtained for the two atomically abrupt defect-free interfaces 
(configurations A and B). Table ^ shows that the p-type SBH is only slightly reduced (by 
0.05 eV, of the same order of magnitude of our numerical precision) in going from the N- to 
the Ga-terminated interface. Similar result JilS were obtained for the GaAs/Al interface, 
where the p-type SBH was found to be reduced by as much as 0.1 eV in going from the As- 
to the Ga-terminated case, in agreement with experiment^l'S. This trend was attributedil 
to the higher electron affinity of the anion-terminated compared to the cation-terminated 
interface, which tends to lower the n-type SBH. 

If we compare the SBH values for the A and C configurations, we note that the p-type 
SBH is almost unaffected by the presence of the AIN intralayer (the 0.05 eV difference being 
of the order of our numerical precision). Let us explain our findings. Since the resulting 
structural relaxations are negligible, the A and C structures differ significantly only for 
the chemical infiuence of the Al atom's presence in the layer closest to the interface. We 
have calculated the difference between the double macroscopic average of the valence charge 
density for the A and C structures (denoted in the following as n^^j'\z)). A large difference in 
the interfacial dipole was not observed due to the similarity of the electronegativity of the Ga 
and Al species. We evaluated its even and odd contribution relative to the substituted cation 
position taken as zero. The even contribution of fi^'j^ {z) represents the chemical difference 
between Ga and Al (e.g. the different atomic wave functions) and, being a quadrupolar 
periodic charge (the monopole term of the electronic part being exactly compensated by the 
nuclear part), that is periodic along x and y but monolayer thick along z, it does not give 
rise to any potential discontinuity which could affect the SBH. On the other hand, the odd 
contribution is expected to represent the charge rearrangement due to chemical substitution, 
which should give rise to charge depletion-accumulation or, equivalently, to an extra-dipole 
that may modify the GaN/Al SBH. Such a dipole can only arise from the asymmetry around 
the substituted cation (due to the presence of the interface), and therefore represents a next 
nearest neighbors effect. In Fig.^, we show the even (solid line) and odd (dashed line) 
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terms of ^^^j// is clear that the odd contribution is very small, relative to the even 

one. Therefore, as expected, Tr^i// i^) ^'^^^ modify the final SBH between GaN and Al. 
(Similar results were obtained from first principles calculations performed on the GaAs/Al 
system0). 

Let us now consider configuration E. Since we have just shown that the presence of an 
interface AIN layer (configuration C) in place of the usual GaN layer (configuration A) does 
not modify the SBH, in a similar way we expect that the presence of an interface Gao.5Alo.5N 
layer will not affect the final SBH as well. This is exactly what happens, as shown by an 
essentially identical SBH for configurations A and E (see Table ^ second row). 

We now discuss the transitivity ruleil for configuration C. In Ref. ^ we showed that 
this simple rule in the abrupt system is fulfilled almost exactly for hypothethically "lat- 
tice matched" systems. In contrast, the difference in lattice constants inherent to the 
GaN/AlN/Al system was seen to produce deviations of about 0.1 eV, due to differences 
in structural interface relaxationJU. Let us now observe that configuration C can be seen as 
a sequence of two interfaces (as far as a single monolayer may represent most of electrostatic 
potential line-up): a GaN/AlN heterojunction and an AIN/Al semiconductor/metal inter- 
face. According to the transitivity rule and recalling the results obtained for the potential dis- 
continuity in the AlN/Al and GaN/AlN casesil, we should therefore obtain: ^Bp{GaN/Al) 
= ^Bp{AlN/Al) - AE^{GaN/AlN) = 1.80 - 0.76 = 1.04 eV, which is in excellent agreement 
with our FLAPW resuh, ^%^{GaN/Al) = 1.07 eV. 

Now consider the D structure; note that its p-type SBH is very similar to the value 
obtained for the other structures. This can be explained considering that the D and B 
structures differ essentially by the isovalent substitution Ga^Al in the semiconductor layer 
closest to the interface. As shown for the A and C structures, the AIN intralayer does not 
alter appreciably the SBH; a similar situation can therefore be expected upon going from 
the B to the D structures, as actually found in Table ||. In more detail, the SBH variation 
on going from B to D has opposite sign, relative to the one seen on going from A to C; this 
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may be related to slightly different atomic relaxations in the two cases. 

Let us now discuss how the SBH can be altered by displacing the atoms in the metallic 
side of the interface. In order to investigate this subject, we make use of the Born effective 
charges (Z*) concept, Z* being the dipole linearly induced by a unitary displacement of a 
single ion in an otherwise perfect crystal. According to Ref. ^ if an atom with Z* is 
displaced by an amount u from its equilibrium position, the resulting SBH change can be 
expressed as: A V = 8 it e'^ Z* u/{eoo a^), where ^oo is the dielectric constant and a is the 
cubic lattice constant of the medium. Of course, it is expected that in a metal, Z* vanishes 
due to perfect screening. In order to evaluate the Z* of the metallic Ga and Al atoms, 
we performed FLAPW calculations for some D type structures, obtained by varying the 
interface distance between the Ga^ and A\} atoms with respect to the equilibrium distance. 
We found, taking u as large as -u = 0.6 a.u., that the variation of the SBH was less than 
0.05 eV. This confirms, as expected, that the Z* for the metallic atoms is nearly vanishing 
and that variations of the SBH as a function of the displacements of the metallic atoms in 
the metal side are therefore negligible. 



IV. CONCLUSIONS 

First principles calculations have been performed for the GaN/Al junction, consider- 
ing various interface morphologies, namely different atomic terminations for GaN, interface 
atomic exchanges and Ga^Ali.^^N intralayers. Our results show that the Ga^Al atomic 
exchange has an enthalpy change of 0.1 eV per atom. 

We find that the first few layers of semiconductor are metallic due to the penetration of 
MIGS. For this reason, the presence of ultra-thin intralayers does not result in a bandgap 
enhancement near the interface. In addition, the effects of the interface configuration on 
the electrical properties of the junction are found to be negligible: the Schottky barrier 
height varies by less than 0.1 eV, changing the interface morphology. In particular, we 
found that the p-type Schottky barrier height is only slightly reduced by the presence of 
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(i) a Ga-termination of the semiconductor side of the junction, (ii) a. Ga<->A1 atomic swap 
at the immediate interface and (Hi) a Gaa;Ali_3.N {x = 0.5, 1) intralayer. Changes in the 
SBH are explained in terms of the different character of the interface bonds in the two Ga- 
and N- terminated atomically abrupt cases and of the Ga2:Ali_j;N (x — 0.5, 1) intralayer 
electronic properties. Our results show good agreement with other theoretical values for the 
p-type SBH obtained for the similar GaAs/Al system. 
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TABLES 



TABLE I. Structural parameters for the configurations examined (all values in atomic units). 



Sys. Interface configuration and interplanar atomic distances 

A Gai 2.12 ^ Nj 2.17 Alj <— 3.67 Al 

B X, ^ 2.2:i Ga, — :i,Oi ^ Al, — :i.(i7 ^ Al 

C Ga^ ^ 2.12 N] ^ 2.12 Al] ^ 2.13 ^ 2.17 Alf ^ 3.67 Al 



D Gai ^ 2.12 N} -S- 2.12 A\} ^ 2.12 ^ Nf 2.12 Gaf <- 3.00 Alf 3.67 ^ Al 

E Gaj ^ 2.12 N] ^ 2.12 ^ Xi ^ 2.12 ^ Nf ^ 2.17 ^ Al^ ^ 3.67 ^ Al 



TABLE II. p-type Schottky barrier heights for the different n GaN + m Al interfaces (values 
in eV). 





A 


B 


C 


D 


E 


n + m 


GaN/Al N-term. 


GaN/Al Ga-term. 


AIN intral. 


Ga<->A1 swap 


Gao.5Alo.5N intral. 


15+7 


1.12 


1.06 


1.07 


1.11 




11+5 


1.30 








1.30 
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FIGURES 

FIG. 1. Different interface geometries considered. Small squares regions indicate the GaN 

substrate; large square regions represent the Al overlayer. Abrupt GaN/Al N-terminated interface 
- system A (panel (a)); abrupt GaN/Al Ga-terminated interface - system B (panel (b)); AIN 
intralayer - system C (panel (c)); system D showing the Ga ^ Al atomic swap (panel (d); GaAlN 
alloy intralayer - system E (panel (e). 

FIG. 2. Valence charge density projected on planes cutting different interface bonds: (a) bulk Al 
atoms in configuration B; (b) Gaj-Alj atoms in configuration B; (c) Ga|-Al| atoms in configuration 
D. Contour lines are spaced by 0.01 electrons/cell. 

FIG. 3. Even (solid line) and odd (dashed line) contributions to n^-jj^(z), defined as the dif- 
ference between the double macroscopic average of the valence charge density for the A and C 
configurations. Units in electrons/basis-area. 



16 



REFERENCES 



^For a review, see W. Monch, Rep. Prog. Phys. 53, 221 (1990). 

^ E. H. Rhoderick and R. H. Williams, Metal-Semiconductor Contacts, Clarendon, Oxford, 
1988. 

^L. J. Brillson, Surf. Sci. Rep. 2, 123 (1983). 

^ J. D. Guo, C. 1. Lin, M. S. Feng, F. M. Pan, G. C. Chi and C. T. Lee, Appl. Phys. Lett. 
68 (2), 235 (1996); Zhifang Fan et al, Appl. Phys. Lett. 68 (12), 1672 (1996). 

^ S. K. Ghandi, VLSI Fabrication Principles: Silicon and Gallium Arsenide, Wiley, New 
York, p. 469, 1983. 

^T. W. Hickmott, P. M. Solomon, R. Fischer and H. Morkoc, J. Appl. Phys. 57, 2844 
(1085); 'W. I. Wang and F. Stern, J. Vac. Sci. Technol. B3, 1280 (1985); J. Batey, S. L. 
Wright and D. J. Dimaria, and T. N. Theis, J. Vac. Sci. Technol. B3, 653 (1985). 

^S. G. Louie, J. R. Chelikowsky and M. L. Cohen, Phys. Rev. B 15, 2154 (1977). 

«A. Zunger, J. Vac. Sci. Technol. 19, 690 (1981); A. Zunger, Phys. Rev. B 24, 4372 (1981). 

9 J. Ihm and J. D. Joannopoulos, J. Vac. Sci. Technol. 21, 340 (1982); J. Ihm and J. D. 
Joannopoulos, Phys. Rev. Lett. 47, 679 (1981); J. Ihm and J. D. Joannopoulos, Phys. 
Rev. Lett. 47, 1568 (1981); J. Ihm and J. D. Joannopoulos, Phys. Rev. B. 24, 4372 (1981). 

° C. A. Swarts, J. J. Barton, W. A. Goddard III and T. C. McGill, J. Vac. Sci. Technol. 17, 

869 (1980); 

^D. J. Chadi and R. Z. Bachrach, J. Vac. Sci. Tech. 16, 1159 (1979). 

2 J. P. A. Charlesworth, R. W. Godby and R. J. Needs, Phys. Rev. Lett. 70, 1685 (1993). 

3 J. Bardi, N. Binggeh and A. Baldereschi, Phys. Rev. B 54, 11102 (1996); C. Berthod, N. 
Binggeh and A. Baldereschi, Europhys. Lett. 36, 67 (1996). 



17 



Mark van Schilfgaarde, E. R. Weber and N. Newman, Phys. Rev. Lett. 73, 581 (1994). 

P. Phatak, N. Newman, P.Dreszer, and E. R. Weber, Phys. Rev. B51, 18003 (1995). 

V. M. Bermudez et a/., J. Appl. Phys. 79, 110 (1996); V. M. Bermudez, J. Appl. Phys. 
80, 1190 (1996). 

P. Skeath, I. Lindau, P. Pianetta, P. W. Chye, C. Y. Su, and W. E. Spicer, J. Electron. 
Spectres.. Relat. Phenom. 17, 259 (1979); P. Skeath, I. Lindau, C. Y. Su and W. E. 
Spicer, Phys. Rev. B 28, 7051 (1983); L. J. Brillson, R. Z. Bachrach, R. S. Bauer and J. 
McMenamin, Phys. Rev. Lett. 42, 397 (1979); R. Z. Bachrach, J. Vac. Sci. Technol. 15, 
1340 (1978). 

T. Kendelewicz, W. G. Petro, LA. Babalola, J. A. Silberman, I. Lindau, and W. E. Spicer, 
J. Vac. Sci. Technol. Bl, 623 (1983). 

S. A. Chambers, Phys. Rev. B 39, 12664 (1989). 

D. C. Sun, H. Sakaki, H. Ohno, Y. Sekiguchi and T. Tanoue, Inst. Phys. Conf. Ser. 63, 
311 (1981); S. P. Svensson, G. Landgren and T. G. Andersson, J. Appl. Phys. 54, 4474 
(1983). 

N. Newman, W. E. Spicer, and E. R. Weber, J. Vac. Sci. Tech. B5, 1020 (1987). 

S. Picozzi, A. Continenza, S. Massidda and A. J. Freeman, Phys. Rev. B 57, 4849 (1998). 

H.J.F.Jansen and A. J. Freeman, Phys. Rev. B 30, 561 (1984); E.Wimmer, H.Krakauer, 
M.Weinert and A. J. Freeman, Phys. Rev. B 24, 864 (1981) and references therein. 

P.Hohenberg and W. Kohn, Phys. Rev. 136, B864 (1984); W.Kohn and L.J.Sham, 
ibid.145, 561 (1966). 

H.J.Monkhorst and J.D.Pack, Phys. Rev. B 13, 5188 (1976). 

Note that, as usual in FLAPW calculations, there is excellent agreement between the cal- 



18 



culated and experimental lattice constant for GaN (a = 8.47 a.u. as reported vaNumerical 
Data and Functional Relationships in Science and Technology, edited by K. H. Hellwege, 
Landolt-Bornstein Tables, Group 111, Vol. 17a (springer. New York, 1982). 

^^A. Di Pomponio, A. Continenza, R. Podloucky and J. Vackar, Phys. Rev. B 53, 9505 
(1996). 

^® Recall that the energy required to form, for example, a (GaN)„/(Al)^ superlattice (SL), 
with n GaN layers and m Al layers starting from the binary constituents can be expressed 
as A if = \Ef^ \ —n-\Ef^^^\ —m - \E^^\ where the total energies of the binaries are obtained 
at the same strain conditions as in the superlattice. 

29 J. Neugebauer, T. Zywietz, M. Scheffler, J. E. Northrup and C. G. Van de Walle, Phys. 
Rev. Lett. 80, 3097 (1998). 

regarding relative 

3°R. G. Dandrea and C. B. Duke, J. Vac. Sci. Techn B 11 (4), 1553 (1993). 

31 T. Ohno, Phys. Rev. B 45, 3516 (1992). 

32 J. S. Foresi and T. D. Moustakas, Appl. Phys. Lett.62, 2859 (1993). 

33 W.I. Wang, J. Vac. Sci. Technol. B 1, 574 /1983). 

34 A. Y. Cho and P. D. Dernier, J. Appl. Phys. 49 3328 (1978). 

35 c. Berthod, N. Binggeh and A. Baldereschi, J. Vac. Sci. Technol. B 14 (4), 3000 (1996). 

3^ G. Margaritondo and P. Perfetti in Hetero junction Band Discontinuities, edited by F. 
Capasso and G. Margaritondo (North- HoUan d, Amsterdam, 1987), p. 59. 

3''' K. Kunc, in Electronic Structure, Dynamics and Quantum Structural Properties of Con- 
densed Matter, edited by J. T. Devreese and P. Van Camp (Plenum, New York, 1985), 
p.227. 



19 



A. Ruini, R. Resta and S. Baroni, Phys. Rev. B 56, 14921 (1997). 



20 



Al 



N 



Ga 



GaN 



N-term. Ga-term. 



AIN 
intral. 



(a) 



(b) 



AIN 




Ga^Al GaAlN 
swap intral. 

(d) (e) 



FIG. 1 - Picozzi et al. 




Fig. 2 - Picozzi et al. 





0.0015 - 




0.0010 - 




0.0005 - 








0.0- 








-0.0005 - 




-0.0010 - 




-0.0015 - 




Ga N Ga N 



X 



N 



Al 



Al 



X = Ga for the A system 
X = Al for the C system 



Fig. 3 - Picozzi et al. 



